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The ground state and low lying excitation of the spin 1/2 Heisenberg ladder with antifer- 
romagnetic leg (J) and ferromagnetic rung (—A J, A > 0) interaction is studied by means of 
the density matrix renormalization group method. It is found that the state remains in the 
Haldane phase even for small A ~ 0.02 suggesting the continuous transition to the gapless 
phase at A = 0. The critical behavior for small A is studied by the finite size scaling analysis. 
The result is consistent with the recent field theoretical prediction. 
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I. INTRODUCTION 



Recently, the spin-1/2 Heisenberg ladder model has been studied extensively 
|l|^,?,plHr7| motivated by its relevance to high superconductivity and Haldane gap 
problem. [1^21] The former is concerned with the fully antiferromagnetic ladder, while the 
latter is related to the ladder with ferromagnetic rung and antiferromagnetic leg interactions. 
In the present paper, we concentrate on the latter model. 

Several years ago, the present author studied this model using the projctor Monte Carlo 
method. [|T| It has been further studied by the non-linear a model, the bosonization 
method ||3|JTT|-[l3| and exact diagonalization. For small interchain coupling, however, 
their conclusions are still controversial. The bosonization calculation by Watanabe, Nomura 
and Takada 0] support the Kosterlitz-Thouless (KT) type behavior with the energy gap 
proportional to exp(— const. /-\/A), and the nonlinear a model analysis predicts the pos- 
siblity of the KT type transition to the gapless state at nonzero positive A. Although the 
projector Monte Carlo simulation |l| seems to support one of these two kinds of possibility, 
the data are rather scattered. On the other hand, the non-Abelian bosonization studies by 
Strong and Millis |^1],|I2| suggest the continuous transition at A = where energy gap is 
proportional to A. Using the same method, Totsuka and Suzuki |]13[ predict the presence of 
the logarithmic correction to the gap as A(ln(l/A))^/^. The exact diagonalization study of 
this model has been also performed by Watanabe. [l^ Nishiyama et al. and Terai 
made the perturbational approach from A = for the finite size system. The density matrix 



renormalization group (DMRG) [|22| , |23| study is also carried out by Narushima et al. fT? 



These works support the continuous transition at A = with the energy gap exponent ^ 1 
with possible logarithmic correction. It is the purpose of the present work to investigete this 
problem further analyzing the DMRG data by the finite size scaling method. 

In the next section, we explain the model and the method of the DMRG calculation. 
The results are analyzed in sections 3 and 4. The last section is devoted to summary and 
discussion. 
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II. HEISENBERG LADDER 



We consider the spin 1/2 Heisenberg ladder consisting of two coupled Heisnberg chains 
represented by the following Hamiltonian H. 

N 

H = 2JY: [Sfst, + SfSf_,, - XSfSf] , (J, A > 0) (2.1) 

i=l 

where S°^{a = A or B) is the spin operator with spin 1/2. The length of the ladder is 
denoted by A^. For A ^ oo, the spins Sf and Sf form a local triplet and this model reduces 
to the spin-1 Heisenberg chain. On the other hand, this model is decoupled into two spin 
1/2 antiferro magnetic Heisenberg chains for A = 0. 

In the DMRG transform, [^,^ the number m of the states retained in the left and right 
half systems ranges from 80 to 140. We use the mixed algorithm of finite size and infinite 
size method. Namely, during the increase of the system size by the infinite size algorithm, 
we sometimes improve the accuracy by moving the boundary between the left and right half 
systems until the eigenvalues converge keeping the system size constant. In order to save 
the computational time, however, we only move the boundary among the sites which are 
added after the last finite size iteration. The data are taken only for the sizes for which the 
finite size iteration is made. For m = 140, only the infinite size algorithm is used unless 
specifically mentioned. On each step of the infinite size iteration, 4 spins are added. In the 
following, we only consider the case of even N. 

In the calculation of the low lying excitation energy eigenvalues, we have constructed 
the density matrices of the left and right half systems from the lowest 3 states of the whole 
system with fixed >S'*°*(= 2;-component of the total spin). Let us denote the energy eigenvalue 
of the n-th excited state with 5**°* by E{n, 5**°*) where n = corresponds to the ground state 
with given value of 5"*°*. For large A^, the calculated energy gap depend on m. We assume 
that the finite m correction behaves as 1/m^ for large m as pointed out by Narushima et 



al. [T^ Taking this into account, we tried two different ways of m-extrapolation. First, the 
data are extrapolated from m = 100, 120 and 140 linearly to 1/m^ by the least squares fit. 
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Secondly, 1/m^-correction is included using the data for m = 80, 100, 120 and 140. The data 
points are the average of these two extrapolated values and the error bars are estimated 
from the difference between them. 



III. KENNEDY TRIPLET 



One of the most remarkable features of the Haldane phase is the presence of the four-fold 
quasi- degeneracy of the ground state and low lying states |^ in the open boundary system. 



In analogy with the S = 1 antiferromagnetic Heisenberg chain, the ground state is a singlet 
and the three lowest excited states form a triplet called Kennedy triplet if the ladder 
length is even. 

Fig. 1 shows the semi-log plot of the size dependence of the excitation energy = 
E{1, 0) — E{0, 0) of the lowest excited states with 5*°* = extrapolated to m ^ oo. The 
convergence of the numerical diagonalization procedure becomes poor for some points due 
to the reason explained in the next section. These data are excluded from the analysis. For 
large A^, the excitation energy vanishes exponentially with the increase of the system size 
within the error bars confirming that the system is in the Haldane phase down to A ~ 0.02. 
This implies that the continuous transition to the gapless phase takes place at A = 0. The 
same conclusion is also obtained by Narushima et al. il^ 



IV. HALDANE GAP 

Taking into account the presence of the Kennedy triplet, the Haldane gap may be es- 
timated by extrapolating the difference /S.Eq = E{2,0) — £"(1,0) |1^ to N oo. The 
system size dependence of AEq is shown in Fig. 2 for A = 0.02,0.05,0.1,0.15 and 0.2 and 
4 < < 102. The second excitation gap AEi = £'(3,0) — £(1,0) is also shown for small 
systems (A^ < 12) which can be diagonalized exactly. From this figure, we find that the gap 
AEq has a maximum at A^ = A^max(A) which increases with the decrease of A. 
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This maximum results from the level crossing which takes place in the finite size system 



at A = Ac(A^). This level crossing is already pointed out by Watanabe [|I^ and the physical 
picture of these excitations are also given. Roughly speaking, we may regard A^max(A) as 
the inverse function of Xc{N), although this is not well-defined in the mathematical sense 
because is a discrete variable. 

The gap AEq for > A^max(A) continue to the second excitation gap AEi for < 
NuiaxW as shown in Fig. 2. Let us define the gap Aq by AEq for > A^max and AEi 
for A^ < A^max, because this state corresponds to the quintuplet excitation. [|l4l In the 
thermodynamic limit A^ — > oo, Aq corresponds to the Haldane gap. Similarly, second triplet 
excitation gap At is defined by AEi for A^ > A^max and AEq for A^ < A^maxj although we 
did not calculate AEi for A^ > A^max because the accuracy of the higher excitation energy 
becomes worse in the DMRG method. Nevertheless, we have shown At for A = 0.05 for large 
A^. These eigenvalues are obtained as the second excited state in the DMRG caluculation. 
However, we interprete these values as At instead of Aq, because it is on the line smoothly 
extrapolated from At. We expect that this occurs because the states important for the 
quintuplet states are already discarded during the iteration for A^ < A'max- Due to the same 
reason, the gap Aq is less accurate than At for large A^. This makes difficult to determine 
the precise value of the exponent u for Aq. 

We therefore start with the analysis of the system size dependence of At by the finite size 
scaling method. In this analysis, we use only the data obtained after the finite size iteration 



even for m = 140. According to Strong and Millis ITTIJT^ and Totsuka and Suzuki, |T3 
the field theoretical calculation using the non-Abelian bosonization technique yields u = 1. 
Totsuka and Suzuki [|13[ also noted the presence of the logarithmic correction to the gap as 



A(ln(l/A))^/^. Expecting the similar behavior for At, let us assume the finite size scaling 
formula for At as, 

A^At(A^) = f{N{\nN)''X'') (4.2) 



This scaling formula is also employed by Terai. ||16l As discussed by Spronken et al. [^, this 



scaling relation implies that the gap At of the infinite system behaves as At ~ A^(ln(l/A))". 
We have tried to fix both a and v by the least squares fit to the universal curve / using 
the data for A = 0.02,0.05,0.1 and 0.2 and 4 < < 102. The data with large error bars 
((5(A^At) >^1) are excluded from the analysis. However, the least squares deviation is rather 
insensitive to the value of a while the optimum value of v remains around unity irrespective 
of the choice of a. For example. Fig. 3 shows the plot with v = 0.88 and a = which 
is the optimum choice with a = 0. Fig. 4 shows the fit with a = 0.4 and u = 1 which 
is the optimum choice with u = 1. In both cases, the data A^At excellently fits the single 
universal curve f{x){x = N{lnN)°'X'^). Although the available value of x for At is limited, 
f{x) behaves linearly for large x within the present data. The scaling plot of Aq is also 
shown. For large A^, it also follows another universal curve which also tends to a linear 
function of x with the same scaling exponent fairly well. This suggests that the critical 
behavior of Aq is also governed by the same exponents as those for At. 

We can obtain almost equally good fit using other sets of values of a and u as far as ~ 1 
both for At and Aq. This clearly excludes the possibility of KT type behavior. However, 
this in turn means that it is difficult to determine the value of a solely from the numerical 
data. Therefore, taking into account the estimations from other works [|TT]-|I^,|T^-[T^ for 
the exponents of the Haldane gap and assuming that Aq and At has the same exponents 
as discussed above, we fix the value of u to unity and search for the optimum value of a. 
In this way, we find a = 0.4 ± 0.07. The error of a is estimated from that of At due to 
m-extrapolation. This value is close to the analytical estimation a = 0.5 by Totsuka and 
Suzuki |jl3| and perturbational estimation of [1^ . 



The finite size scaling plot in Fig. 4 also suggests that the level crossing point is also 
scale invariant as Ac(A^)A^(lnA^)" = const. ~ 5.2. This implies that Ac tends to as A^ oo, 
confirming that the gap Aq remains the lowest gap down to A = in the thermodynamic 
limit. 
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V. SUMMARY AND DISCUSSION. 



We have studied the low energy excitation of the spin-1/2 Heisenberg ladder with fer- 
romagnetic rungs and antiferromagnetic legs using the DMRG method. It is found that 
the ground state remains in the Haldane phase down to A = 0.02 suggesting the continous 
transition to the gap less phase at A = 0. The exponent v of the second triplet gap At is 
estimated to be around unity using the finite size scaling analysis. Fixing z/ = 1, we find the 
logarithmic correction as At ~ A(ln(l/A))" with a = 0.4 ± 0.07. Although the numerically 
obtained value of Aq is not enough accurate to determine the exponents independently from 
At, it is expected Aq also shows the same critical behavior taking into account the fact 
that Aq also fits to the universal curve fairly well with the same scaling within the error 
bars. It is also shown that Aq remains the lowest gap (Haldane gap) down to A ^ in the 
thermodynamic limit. 

Several years ago, the present author proposed the possibility of the transition at finite 
A based on the projector Monte Carlo simulation and the mapping onto the nonlin- 
ear a model. However, the Monte Carlo data are rather scattered and the nonlinear a 
model analysis is based on the semiclassical approximation and the artificially introduced 
anisotropy. Therefore, we believe the present result is much more reliable. Actually, the 
non-Abelian bosonization analysis employed by Strong and Millis |n|JI31 Totsuka and 



Suzuki, []13[ which leads to the results consistent with our calculation, may be regarded as 
the refined version of the nonlinear a model analysis suited for the isotropic case. 
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FIGURES 

FIG. 1. The system size dependence of the energy gap /S.k between the ground state and the 
Kennedy triplet for A = 0.02,0.05,0.1,0.15 and 0.2. 

FIG. 2. The system size dependence of the energy gap AE'o (open symbols) and AE^i (filled 
symbols) for A = 0.02,0.05,0.1,0.15 and 0.2. 

FIG. 3. The finite size scaling plot with u = 0.88 of the energy gap At (open symbols) and Aq 
(filled symbols) for A = 0.02,0.05,0.1,0.15 and 0.2. 

FIG. 4. The finite size scaling plot with v = 1 and a = 0.4 of the energy gap At (open symbols) 
and Aq (filled symbols) for A = 0.02,0.05,0.1,0.15 and 0.2. 
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